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Abstract 

We show that the spectrum of the three dimensional c/> 4 theory 
in the broken symmetry phase contains non-perturbative states. We 
determine the spectrum using a new variational technique based on 
the introduction of operators corresponding to different length scales. 
The presence of non-perturbative states accounts for the discrepancy 
between Monte Carlo and perturbative results for the universal ra- 
tio £,/£,2nd- We introduce and study some universal amplitude ratios 
related to the overlap of the spin operator with the states of the spec- 
trum. The analysis is performed for the </> 4 theory regularized on a lat- 
tice and for the Ising model. This is a nice verification of the fact that 
universality reaches far beyond critical exponents. Finally, we show 
that the spectrum of the model, including non-perturbative states, ac- 
curately matches the glueball spectrum in the 71*2 gauge model, which 
is related to the Ising model through a duality transformation. 
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1 Introduction 



The idea of universality and the perturbative analysis of 4 theories have 
proved to be very powerful tools in the study of several statistical mechan- 
ics models (see e.g. Ref. [||] and references therein). In particular, with 
the advent of the last generation of high precision simulations for the three 
dimensional Ising model, an impressive agreement has been found between 
numerical results and perturbative predictions from 4 theory. Almost all 
universal quantities, both critical indices and amplitude ratios, agree within 
error bars ||. With one small exception. The ratio between the "true" 
exponential correlation length £ and the second moment one ^nd which is 
predicted to be £/£2nd ~ 1.0065, from a <p A calculation in the broken symmet- 
ric phase H |J], turns out to have a larger value C/Cnd ~ 1.031(6) both from 
Montecarlo simulations and from low temperature expansions in the three di- 
mensional Ising model. This paper deals with this discrepancy. We shall show 
that it can be understood as due to the presence of new non-perturbative 
states in the spectrum of the theory. We shall also evaluate various universal 
ratios involving masses and overlap constants of the non-perturbative states 
of the spectrum. 

The analysis of the spectrum has been performed by a variational method, 
by introducing a set of operators analogous to the one introduced in Ref. || 
in the context of lattice gauge theories. These operators correspond to dif- 
ferent length scales and are constructed recursively. We were able to identify 
precisely the first states of the spectrum. This new method (which is not 
restricted to the Ising case, but can be easily extended to any spin model) is 
among the main results of this paper and we shall discuss it in great detail. 

We shall then show that universality holds also for this non trivial part 
of the spectrum, by directly simulating the lattice version of the 4 theory 
and again finding the same pattern of non-perturbative states and the same 
values of the universal ratios. Finally, we shall show, by using duality in the 
Ising model that these new states coincide with the lowest excitations of the 
glueball spectrum of the (dual) 7L 2 gauge model. 

We will also present a detailed study of some universal amplitude ratios 
related to the overlap of the spin operator with the low-lying states of the 
transfer matrix. In particular, we will define and study a universal amplitude 
R which involves the overlap of the spin operator on the lowest state: for this 
quantity we can compare numerical results with perturbative calculations, 
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which we extend to two loop level. 

This paper is organized as follows: in Sec. 2 we introduce the Ising and 4> 4 
models and the observables we will be interested in. In Sec. 3 and Sec. 4 we 
introduce two universal quantities ^/^nd and R (related to the overlap con- 
stant of the spin operator) and discuss the existing numerical and analytical 
results about these two quantities, including their perturbative evaluation 
which we have extended to two loop level. In Sec. 5 we describe the new 
variational method we have used to determine the spectrum. Sec. 6 contains 
our Monte Carlo results for the spectrum of the model obtained with the 
variational approach: these results show unambiguously the existence of a 
non-perturbative state. In Sec. 7 we discuss duality and the relationship 
with the glueball spectrum of the 7L 2 gauge model. In Sec. 8 the spin- 
spin correlation function is reconsidered taking into account the existence 
of non-perturbative states in the spectrum; a new universal quantity R 2 is 
introduced, which is related to the overlap constant of the spin operator on 
the lowest non-perturbative state. Finally Sec. 9 is devoted to some conclud- 
ing remarks, and in the Appendix we collect the details of the perturbative 
calculations. 

2 The models and the observables 

2.1 Ising model 

The Ising model is defined by the action 



where the field variable s n takes the values —1 and +1; n = (no, m, n 2 ) 
labels the sites of a simple cubic lattice and the notation < n,m > indicates 
that the sum is taken on pairs of nearest neighbour sites only. The coupling (3 
is proportional to the inverse temperature, (3 = t^. We shall consider in the 
following m and n 2 as "space" directions and no as the "time" direction and 
shall sometimes denote the time coordinate no with r. The high and low T 
phases are separated by a critical point at a coupling whose value is known 
with very high precision M: (3 C = 0.2216543(2)(2). A peculiar property 
of the Ising model, which will play an important role in the following, is 




(1) 



<n,m> 
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the existence of an exact duality transformation which relates it to the H12 
gauge model. This transformation is known as Kramers-Wannier duality 
and relates the two partition functions: 

£ = -ilog[tanh(/3)] , (2) 

where (3 will be denoted in the following as the "dual coupling". Using the 
duality transformation it is possible to build a one-to-one mapping of phys- 
ical observables of the gauge system into the corresponding spin quantities. 
In particular the inverse of the mass of the lowest state in the spectrum of 
the gauge model ^ gaU ge is mapped into the correlation length of the spin Ising 
model £. 

2.2 4 model 

The lattice version of the </> 4 model is given by the action 

<n,m> n n 

where now the field variable n assumes all possible real values. In the limit 
A — > 00 the standard Ising model is recovered. In the space spanned by 
the two coupling constants (/3, A) the model has a second order critical line 
which belongs to the same universality class as the Ising model. In fact the 
two models share the same 7L2 global symmetry which is broken in the low 
T phase. A peculiar feature of this model is that by suitably tuning A one 
can reach a point in which the corrections to scaling (proportional to £~ w ) 
disappear (see Refs. || |7j). It turns out that the optimal value is A = 1.1. 
With this choice the critical coupling is (3 = 0.3750966(4). 

We shall be interested in the exponential and second moment correlation 
lengths in the broken symmetry phase. Their definition is the same for the 
Ising and 4 models. 
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2.3 Exponential correlation length. 

The exponential correlation length is defined in terms of the long distance 
behavior of the connected two point function: 

1 1 

- = - Jim — log(s d s H )c • (4) 

£ \n\— >oo 1 71 1 

Note that for the 4 model s has to be replaced by 0. It is convenient to 
study the so called time slice correlation functions: The magnetization of a 
time slice is given by 

^no = J~2 ^ j s (n ,ni,n2) ■ (5) 
ni,n 2 

Let us define the time slice correlation function 

G(r) = J2{(Sn S n()+T }-(S no } 2 } . (6) 



no 



The large distance behavior of G(r) is given by 

G(t) = c eM-r/O , (7) 

where £ is the exponential correlation length, and c is a constant. The 
constant c is related to the overlap of the spin operator with the eigenstate 
of the transfer matrix corresponding to the lowest eigenvalue. A universal 
ratio involving the constant c will be discussed in Sec. 4. 

It is important to stress that for small values of r one expects deviations 
from this asymptotic behavior. These corrections are due to higher masses in 
the spectrum. A common tool to extract £ is the so called effective correlation 
length 

ee//(r) = ln(G(r + l))-ln(G'(r)) ' (8) 

The effective correlation length is a monotonically increasing function of r. 
In the limit r — > oo it converges to the exponential correlation length £. 
Notice that in a Monte Carlo simulation with increasing r also the statistical 
errors of £ e // become larger. Therefore a fast convergence is important for 
the numerical study. It turns out that £ e // in the broken phase of the Ising 
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Fig. 1: Data for^ eff (r) at (3 = 0.2275 (taken from Ref. gj). The dashed 
line corresponds to the asymptotic value £ = 2.592 (see Tab. 1). The 
dotted line corresponds to the function £, per t defined in Sec. 3.4 

model is not a good estimator of £, since it requires distances larger than 3£ 
to approach £ with a relative accuracy of 1% (see Ref. |2j for a discussion of 
this point). 

This is clearly visible in Figs. 1 and 2 where data taken from || are 
plotted. In particular in Fig. 1 we have shown the data for £ e ff corresponding 
to (3 = 0.2275 while in Fig. 2 all the data of Ref. are plotted together after 
a suitable rescaling. It is possible to see from Fig. 2 that all the data show the 
same (non-asymptotic) behavior in the range £ < r < 3£. This shows that 
scaling is fulfilled in this range and that the deviation from the asymptotic 
behavior is certainly due to some physical reason (namely to the presence of 
nearby masses in the spectrum) and not to lattice artifacts. 
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Fig. 2: £ e //0") for (3 = 0.23910,0.23142,0.2275,0.2260,0.2240. All 
the data are taken from Ref. j^/. Both ^ e jf and the distance r are 
normalized, for each (3, to the asymptotic value £. 
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2.4 Second moment correlation length. 



The square of the second moment correlation length is defined for a d- 
dimensional model by 

where 

Ho = lim — (s m s n ) c (10) 

L^oo V — ' 
m,n 

and 

/i 2 = lim - V(m- n) 2 (s m s n ) c , (11) 

m,n 

where the connected part of the correlation function is defined by 

(SmSn)c {SmSn) (^m) (12) 

and V is the lattice volume. 

This estimator for the correlation length is very popular since its numer- 
ical evaluation (say in Monte Carlo simulations) is simpler than that of the 
exponential correlation length. Moreover it is the length scale which is di- 
rectly observed in scattering experiments. However it is important to stress 
that it is not exactly equivalent to the exponential correlation length. We 
shall discuss the relation between the two in the next section. 



3 The f/f 2 nd ratio 
3.1 £ versus &nd 

The relation between £ and ^nd can be obtained by noticing that we can 
rewrite /i 2 as follows 

n;m 
^ d-l 

= ^ ~ ( S ™ S ")c ■ ( 13 ) 

n;m 
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Due to the exponential decay of the correlation function this sum is conver- 
gent and we can commute the spatial summation with the summation over 
configurations so as to obtain 



H 2 = dJ2 r * (SO Sr)c , (14) 
r= — oo 

with S no given by Eq. (|). Analogously one obtains 

oo 

fl = ^ ( S S r)c ■ (15) 
r= — oo 

If we now insert these results in Eq. (^), we obtain 

2 2^r=-oo G ( T ) 

Assuming a multiple exponential decay for G(t), 

(S S T ) C oc ^ a exp(-|r|/^) , (17) 

i 

and replacing the summation by an integration over r we get 

£2 = 1 J T =o dr r2 E»c t exp(-r/&) = £ . c^ 3 ^ 

which is equal to £ 2 if only one state contributes. 

An interesting consequence of this analysis is that the deviation of the 
ratio £ / ^nd from the value 1 gives an idea of the density of the lowest states 
of the spectrum. Note that Q > (see Eq. (pE0|) ) and therefore £/^2nd > 1- 
If these are well separated the ratio will be very close to one, while a ratio 
significantly higher than one will indicate a denser distribution of states. 

Notice that besides a discrete sum of exponentials (each corresponding 
to a pole in the Fourier transform of G(r)) we can also have an integral over 
a continuous set of exponential functions, which corresponds to a cut in the 
Fourier transform. This is the case, for instance, for the 4 model above 
the pair production threshold at p — 2m. These cuts can be thought of as 
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the coalescence of infinitely nearby exponentials, and actually, on a finite 
lattice, this is their correct description, since the transfer matrix has only a 
finite number of eigenvalues. The effect of these cuts is, as it happens for the 
isolated states, to enhance the ratio £/£2«d- We shall refer to the contribution 
to the connected correlator due to these terms as the "cut contribution" . 

3.2 MC estimates 

Rather precise estimates for £ 2 nd can be found in Ref. Q for various values of 
j3 in the scaling region of the 3D Ising model. In the same paper estimates 
for £ also appear. However, while ^2nd can be extracted rather easily from 
MC simulations, the determination of £ is much more delicate since, as we 
discussed above, it requires the identification of an asymptotic exponential 
decay. If other states besides the lowest one are present in the spectrum 
of the theory (and this is exactly the situation in which we are interested 
here) they can shadow the asymptotic behavior that we are looking for and 
produce systematic errors in the estimate of £ (this problem was discussed 
in great detail in Ref. @). 

A natural way out of this problem is to use a variational analysis to 
separate the various states in the spectrum (see 0). However this requires 
a wide set of operators with a good overlap with the low-lying states of the 
theory. Constructing these operators is a highly non trivial task in the Ising 
model and has never been attempted up to now (we shall address and solve 
this problem in Sec. 5 below). Fortunately in the case of the 3D Ising model 
there is a nice way to avoid this obstacle. By using duality we can map the 
Ising model into an equivalent H2 gauge model. In the gauge model one 
can easily construct, by looking at Wilson loops (namely products of gauge 
variables along closed contours) of various size and shape, a set of operators 
fulfilling the above requirements and perform a precise variational analysis 
of the spectrum. This route was followed in Ref. M leading to the results 
reported in Tab. 1. From them we extract the estimate £/^2nd = 1.031(6) 
obtained ignoring the first value, which is too far from the critical region, 
and taking a weighted sum of the remaining values (hence assuming that the 
correction to scaling terms are negligible within the errors for this ratio). 



9 



Table 1: Comparison of the results for the exponential correlation length 
obtained directly in the Ising model without a variational analysis (fourth 
column, data taken from Ref fflj, denoted as £ e ff) an d those obtained in the 
dual 7L 2 gauge model with a variational analysis (third column, data taken 
from Ref denoted as i ga uge)- P denotes the dual of f3 and is reported for 
completeness in the first column. In the hfth column we report the values for 
£,2nd obtained in Ref. fflj and hnally in the last column the ratio £, gaU ge/£,2nd 
that we consider as our best estimate for ^/^ 2n d- 



p 


P 


^gauge 




^2nd 


£,gauge/ ^2nd 


0.72484 


0.23910 


1.296(3) 


1.2851(28) 


1.2335(15) 


1.051(4) 


0.74057 


0.23142 


1.864(5) 


1.8637(45) 


1.8045(21) 


1.033(4) 


0.74883 


0.22750 


2.592(5) 


2.578(7) 


2.5114(31) 


1.032(3) 


0.75202 


0.22600 


3.135(9) 


3.103(7) 


3.0340(32) 


1.033(4) 


0.75632 


0.22400 


4.64(3) 


4.606(13) 


4.509(6) 


1.029(8) 



3.3 Series expansion 

It is possible to obtain an estimate for the ratio ^/^nd by using the low 
temperature series published in Ref. [[J. This analysis has been recently 
performed in Ref. ]l(J with the result ^f&nd — 1-031(5) . 



3.4 Perturbative result 

In the framework of the 4 theory we naively expect only one mass in the 
spectrum (hence only one isolated exponential in G(r)). However it is easy 
to see that this result is true only at tree level and that at one loop, a cut 
appears in the Fourier transform of the propagator, starting from twice the 
value of the fundamental mass. As mentioned above this exactly coincides 
with the pair production threshold. The corresponding expression for G(r) 
(which we shall call in the following G(r)p er t) is M: 



G{r) 



pert 



1 

2m ph L 2 



-m ph i 



1 + 



32 4vr 



10 



3u R f°° , e 



16nL 2 m p h 



-/AT 



d\i 2 • (19) 



2m ph fl* - 



III- , 



where u R denotes the dimensionless renormalized coupling, m p h is the phys- 
ical mass, defined as the location of the zero of the inverse correlator in 
momentum space G~ 1 (p) and coinciding with the inverse of the exponential 
correlation length £ (for all the details and definitions concerning perturba- 
tive results, see the Appendix). 

Eq. (ID shows that even perturbatively the long distance behavior of 
time slice correlations is not purely exponential: this, as discussed in the 
previous section, implies a ratio £/^2nd different from 1. Indeed, defining 
the renormalized mass m R as the inverse propagator in momentum space at 
p = 0, so that m R = l/^,2nd, one finds at one loop |3], || 

£ m R u R ( 13 3 



. (20) 



The fixed point value u* R of the renormalized coupling has been recently 
estimated in a high precision Monte Carlo simulation || to be u* R = 14.3(1). 
Plugging this result into Eq. (|20|) we obtain the universal value 



" 1.00652(3) , (21) 



a result which is rather far from the Monte Carlo one. 

A similar exercise is to define the function £ pert by inserting G pert (r) in 
Eq. (^|) and compare it to the Monte Carlo results for £ e //. This comparison 
is shown in Fig. 1. Again, it appears that the perturbative contribution alone 
is not enough to justify the deviation of £ e ff from its asymptotic value. 

We have extended the calculation of £/^2nd to two loop level; details are 
given in the Appendix. The result is 



1 + 0.00573— + 0.00474— . (22) 



6nd " 4tt ' " 16tt 2 

If we plug in the fixed point value u* R = 14.3 we obtain 



> 1.01266 . (23) 

Z2nd 
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One can see that the convergence properties of the series are not very satis- 
factory. Since at the fixed point ur/Ati = 1.14, the two loop contribution is 
as big as 94% of the one loop term. Therefore this result must be taken with 
great caution. Taking it at face value, we see that the two loop contribution 
goes in the right direction, but is not sufficient to close the gap between 
perturbative and Monte Carlo results. 

These results seem to indicate that perturbative effects alone cannot ex- 
plain the rather large value of £/£2n<2 found in numerical simulations or, equiv- 
alently, the preasymptotic behavior of £e//(i")- An unambiguous indication 
that such non-perturbative physics is actually present will be given by the 
determination of the spectrum of the transfer matrix by a variational method 
in Sec. 6. 



4 Universal ratios of overlap amplitudes. 

From the large distance behavior of the correlation function G(t) we can ex- 
tract, besides the ti/^nd ratio, another universal quantity defined as follows. 
If we parameterize the large distance behavior of G{t) as 

G(t) = c exp(-r/0 , (24) 

then the ratio 

(where L is the size of the lattice in the spacelike directions and M denotes 
the magnetization) is universal. 

Since this amplitude combination is not among those usually discussed 
in the literature it is worthwhile to describe its meaning in more detail in 
the framework of quantum field theory and its relationship with the so called 
"overlap constants" . 

Let us look, as an example, at 4 theory in three dimensions at tree 
level in the broken symmetry phase. The point — point connected correlator 
(normalized to the square of the mean magnetization) in momentum space 
is given by 

Wjgc = W f d 3 p j**-v) 

(0(O)) 2 m J (27r)V + m 2 1 ' 
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where denotes the projection of the field on the momentum basis and 
hence is usually referred to as the overlap of the field with the particle 
state of mass m. 

Let us now turn to the slice-slice correlator. The slice operator is defined 

as 

S(t) = -j^ J dx 1 dx 2 (p(x 1 ,X2,t) . (27) 

Plugging this definition into Eq. (|26|) and performing the various integrals 
we find for connected slice-slice correlator G(t): 

G(t) = (S(0)S(r)) c = (S(0)} 2 ^- 2 e- m ^ . (28) 



Hence if we assume the large distance behavior of Eq. we immediately 
recognize that the amplitude ratio defined in Eq. (|25|) exactly coincides with 

This analysis can be straightforwardly extended to the case in which more 
than one massive state is present in the spectrum of the theory. If we expect 
(as in Eq. ( |T7D ) a multiple exponential decay for G(t), 

(S S T ) C oc >^ a exp(-|r|/&) , (29) 

then we can define an independent universal ratio for each state of the spec- 
trum 

* = (30 » 

which coincides with the (square of the) overlap amplitude Fl, that is with 
the projections of the <fi operator on the i th massive state of the spectrum. 

The F^'s encode much interesting information on the theory. In some 
cases, for instance for two dimensional integrable models, they can be eval- 
uated exactly from the S-matrix of the models. 

In the case of the three dimensional Ising model, no exact result is known, 
but it is possible to extract a perturbative estimate of R using the same 
results discussed in Sec. 3.4. We shall deal with this calculation in the next 
section. 
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4.1 Perturbative result 



Looking at Eq. (|T9| ) we see that the constant c in front of the exponential is 
given at one loop by 



1 



2m ph L 2 . 



1 + 



ur 

1287T 



(31) 



Plugging this result into Eq. 
coupling (see Appendix) 



Ur 



using the definition of the renormalized 

M 2 ' 



we obtain 



6 



1 - 



Ur 

1287T 



so that at the fixed point u* R = 14.3 

R ~ 2.298 



(32) 
(33) 

(34) 



Also this calculation can be extended to two loop level as shown in the 
Appendix; as in the case of the amplitude ratios ^/^nd we find that the 
perturbative series has poor convergence properties. The result is 



R 



— ( 1 - 0.03125— - 0.02567^ 
6 



tr, 



4vr 16vr 2 
so that plugging in the fixed point value u* R = 14.3 we obtain 

2.219 . 



(35) 



(36) 



4.2 MC estimate 

It is easy to construct an estimator for the ratio R if we assume the value 
of the correlation length as an input. From Eqs. (0,^) we see that the 
quantity 

(L\ 2 G{r) exp(r/Q _ 
Reff = [j ) W (37) 

converges to R for r — > oo. The behavior of R e ff as a function of r is similar 
to that of £ e //. At short distances it is a decreasing function of r and then 
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reaches a stable plateau for r > 3£. This is clearly visible in Fig. 3 where data 
obtained for four values of (3 (the data are taken from Ref. ||) are plotted 
together. It is interesting to notice that in the region £ < r < 3£ all of the 
four samples follow the same curve, showing that the preasymptotic behavior 
of R is, like the one of £ e //, a physical effect rather than a lattice artifact. 
As in the case of £ e // ; this behavior signals that the correlation function is 
not a single decaying exponential, hence the presence of higher states in the 
spectrum and interaction effects (cuts in the Fourier transform). We shall 
come back to this point in Sec. 7 below. We report in Tab. 2 the values of 
R(t = 3£) which we consider as our best estimates for the asymptotic value 
of R. Let us briefly comment on how we obtained the numbers plotted in 



Table 2: The ratio R for various values of (3. 



p 


R 


0.23142 
0.22750 
0.22600 
0.22400 


2.072(24) 
2.068(18) 
2.060(30) 
2.082(50) 



Fig. 3 and reported in Tab. 2. We used as input data the estimates of £ 
reported in the third column of Tab. 1 (extracted from Ref. ||) and used for 
G(t) the high precision values obtained in 0. The main source of error in 
R e ff comes from £. Due to the exponential in Eq. ([37]) it increases as a linear 
function of r. This explains the rather large errors quoted in Tab. 2. As can 
be seen in Tab. 2 the estimates of R are stable within the errors as a function 
of f3. A naive extrapolation (neglecting corrections to scaling) suggests the 
value R = 2.07(4) at the fixed point, which shows a 10% and 7% deviation 
from the one and two loop perturbative estimates respectively. 

As discussed in Sec. 3, the two loop perturbative evaluation of R is not 
very reliable due to its poor convergence. If one chooses to take it seri- 
ously, one is left with another discrepancy between two loop perturbative 
predictions and numerical results, that is another hint to the existence of 
non-perturbative states in the spectrum. 
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Fig. 3: Reff{r) for (3 = 0.23142 (diamonds), 0.2275 (pluses), 0.2260 
(squares) and 0.2240 (crosses). All the data are taken from Ref. ^ . 
The distance t is normalized, for each (3, to the asymptotic value To 
avoid confusion only the error bars for the (3 = 0.2275 data are reported, 
the other error bars are of the same size. The dashed line corresponds 
to the asymptotic value R = 2.07 quoted in Sec. 4.2 
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5 The variational method for the determina- 
tion of the spectrum 

As mentioned above, the only way to obtain reliable values for the masses 
of a complex spectrum like the one in the broken phase of the Ising model 
is to use a variational technique. Since this is one of the main points of 
our analysis we shall devote this section to a detailed discussion, first of the 
general features of the approach and second of the choice of the operators. 



5.1 Variational analysis 

In general the mass spectrum of a theory is given by the eigenvalues of the 
Hamiltonian H. On the lattice, one diagonalizes the transfer matrix T, which 
is the discrete version of e~ H . For a finite lattice the transfer matrix of the 
Ising model is a real symmetric matrix. Therefore it can be diagonalized. 
Let us denote the resulting eigenvalues by Aj. Then the mass-spectrum is 
given by 

mi = -log(^) , (38) 

where Ao is the largest eigenvalue of T. The basic strategy to evaluate these 
eigenvalues is to compute expectation values of certain correlation functions. 
Masses can then be determined from the decay of these correlation functions 
with the separation in time. 

n t \ /Amw, w (0\AT T B\0) 
G AB (r) = (A(0)B(t)) = | o|rT|o | 

~ (0\A\i)(i\TT\j){j\B\0) =J2 c? B (^ 

= ^cf B exp(-m l r) , (39) 

i 

where \i) denotes the eigenstates of the transfer matrix and[| 

cf_ = (0|A|z)(z|S|0) . (40) 

1 If A = B = S(t) (where S(t) is the slice operators defined above) then the constants 
cf B coincide with the overalp constants a defined in sect. 4 . 
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The main problem in the numerical determination of masses is to find 
operators A and B that have a good overlap with a single state i.e such 
that Cj is large compared with cj, j 7^ i. A first simplification is obtained by 
using the so called "zero momentum" operators, namely operators obtained 
by summing over a slice orthogonal to the time direction, so that all Cj's 
that correspond to nonvanishing momentum vanish. The zero momentum 
operators are just the time slice averages introduced in Sec. 2. 

A systematic way to further improve the overlap is to simultaneously 
study the correlators among several operators A a . One must then mea- 
sure all of the correlations among these operators and construct the cross- 
correlation matrix defined as: 

C a ^T) = (A a (r)A^0))-(A a (r))(A^0)) . (41) 

By diagonalizing the cross-correlation matrix one can then obtain the mass 
spectrum. 

This method can be further improved |J by studying the generalized 
eigenvalue problem 

C(t)4> = \(t,t )C(t )iP , (42) 

where r is small and fixed (say, r = 0). Then it can be shown that the 
various masses m 8 are related to the generalized eigenvalues as follows |J: 

where both r and r should be chosen as large as possible, r >> t and as 
r is varied the value of m; must be stable within the errors. Practically we 
are forced to keep r = to avoid too large statistical fluctuations and at 
the same time r is generally forced to stay in the range r — 1 to 5. This 
method is clearly discussed in Refs. || , to which we refer for further details. 
All the results that we shall list in the next section have been obtained with 
this improved method. In order to give some information on the reliability 
of the estimates we shall also list, besides the numerical values of the masses, 
the value of r at which they have been evaluated. 

5.2 The Operator Basis 

While the formalism of the variational approach to compute the spectrum is 
quite general, a suitable set of operators to compute the correlation functions 
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has to be found for each model separately. Suitable in this context means that 
the wavefunctions of the small mass states are given to a good approximation 
by some linear combination of the operators selected. 

In the case of the 7L 2 gauge model a good set operators for the + channel 
is formed by Wilson loops of various sizes |§ . We kept in the present analysis 
the basic idea that different operators should correspond to different length 
scales, and we included the standard time-slice magnetization in our basis of 
operators. 

We came up with a recursive definition of the operators. The start- 
ing point is the field </>n°o,n 1 ,n 2 — 0n o ,ni,n 2 as h is generated by the Monte 
Carlo. In the case of the Ising model we similarly start with the defini- 
tion 0n°o,ni,n 2 = s n0inijn2 . Then the field (for each time-slice separately) is 
transformed according to the following rule: 

^notnL = sign(«) ((l-w)\u\+wy)) (44) 

with 

v=r An) , /-, \ ±(An) ,A n ) ,An) ,jM x 

T y no,ni,n2 1 v / ^ vr no,ni— l,ri2 rnojfii+ljWa ^'^o^ n i: n 2 — 1 ' no,rii,n2+l/ 

(45) 

w, y and r being free parameters of the transformation. The operators are 
then given by the sum of the <p n over a given time-slice (zero momentum 
projection) 

SW(no) = 5>S«,™ ■ (46) 

ni,n 2 

The correlation matrix is then built by 

Cy(r) =< S(n )®S{n +r)W >-< \S(n )^\ >< \S(n +r)^\ > . (47) 

We performed test simulations with different choices for the parameters w, y 
and r for small correlation length both in the Ising and 4> 4 models. The qual- 
ity of the resulting operator basis was judged by looking at the convergence 
of effective masses towards their asymptotic values. It turned out that r = 
is a good choice for 4 (for Ising we have to choose r slightly different from 
to have a well defined transformation). Also there is no sharp constraint 
on the value of y. Mostly we have taken y equal to the magnetization. In 
the case of w it turned out that values close to are a good choice. Note 
however that taking w exactly equal to would result in all being equal. 
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6 The spectrum: Monte Carlo results 



We simulated the Ising model in its low temperature phase at (3 = 0.23142 
and 0.2275, and the 4 model at fixed A = 1.1 for the three values j3 = 0.405, 
0.385 and 0.3798. The two critical values of (3 are (3 C = 0.2216543(2) (2) 
and p c = 0.3750966(4) for the Ising model and the 4 theory at A = 1.1 
respectively. In both cases we used a single cluster algorithm |lTfl . For the 
application of cluster algorithms in 4 type models see Ref. [12]. 

We used cubic lattices with periodic boundary conditions and size L in 
both spatial directions (those which define the plane in which we perform 
the zero momentum projection of our observables) and chose a size 2L in the 
"time" direction in which we evaluate correlators. L was always chosen such 
that L > 16£. Some information on the simulations is collected in Tab. 2. 

For the whole set of (zero momentum projected) improved operators 
S^ n \n ) we computed the correlators Cjj(r) for all possible translations in 
the time direction. For the values of (3 closer to the critical point, we used 
the set of operators obtained by iterating twice the smoothing procedure de- 
scribed in the previous section. In all five simulations we found the same 



Table 3: Some information on the run in the low temperature phase of the 
Ising model. N op denotes the number of operators used in the variational 
analysis, N srnoothing the number of smoothing iterations between two opera- 
tors. 



p 


L 


measures 


sweeps / measure 


smoothing 


N 

ly op 


0.23142 


30 2 x 60 


300000 


20 


1 


20 


0.2275 


45 2 x 90 


500000 


30 


2 


20 



pattern: 

a] The data for the lowest mass are of much better quality than those 
extracted from the simple spin-spin correlator. This means in particular 
that we have a very precocious approach to the asymptotic result and 
always find a stable plateau starting from values of r of the order of (or 
even smaller than) the correlation length. The quality of the data can 
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Table 4: Some information on the run in the low temperature phase of the 
4 model at A = 1.1. N op and N smoothing are defined as above. 



p 


L 


measures 


sweeps / measure 


-^smoothing 


N op 


0.405 


20 2 x 40 


2300000 


10 


1 


12 


0.385 


40 2 x 80 


850000 


15 


2 


20 


0.3798 


60 2 x 120 


250000 


20 


2 


20 



be better appreciated by looking at Fig. 4 where our data and those 
extracted from the pure spin-spin correlator are compared. 

b] In all cases we can detect the first excitation above the lowest mass and 

evaluate its mass with good precision. For the simulations farther from 
the critical point we can even detect the next excitation (see however 
the next section for a comment on the interpretation of this further 

state) 

c] In applying the variational approach discussed above we always chose 

ro = while the distance r at which the mass is measured varies from 
sample to sample and is reported for completeness in Tabs. 5 and 6. 

d] If we decrease the number of operators involved in the variational anal- 

ysis our data smoothly reach those extracted from the pure spin-spin 
correlator. The rate of this approach gives an idea of the number of 
operators (and smoothing iterations) needed to obtain a good overlap 
with the states of the spectrum. 

e] On finite lattices there is no spontaneous symmetry breaking, because of 

finite action tunneling solutions. The effect of these solutions on the 
mass spectrum is to split each level into a nearly degenerate doublet 
(see e.g. Ref. |13| for a discussion of this point). This splitting is of 



order y/ae~ aL , where a is the interface tension and L the transverse 
size of the lattice. With our choices of parameters we have <jL 2 ~ 30 
(cfr. Ref. [[L4j] where the relevant values of a are reported), therefore 



the splitting cannot be observed with our resolution. 
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2.65 



2.6 



2-55 - 

*eff 2 - 5 " 
2.45 - 

2.4 - 

2.35 ' 1 1 1 1 1 1 

0.5 1 1.5 2 2.5 3 3.5 

Fig. 4: Values of£ e ff obtained with the variational method (pluses) and 
with the standard spin-spin correlator (diamonds) for (3 = 0.2275. The 
dashed line denotes the asymptotic value £ = 2.592 (see Tab. l).The 
distance r is measured in units of £. 
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The values for the (inverse of the) masses that we found are reported in 
Tabs. 5 and 6. In order to address the issues of scaling and universality we 



Table 5: Correlation lengths extracted with the variational approach in the 
3d Ising model. Below each mass we report the value of r at which it has been 
evaluated, tq = is always assumed. The question marks denote the fact 
that the corresponding states are not yet stable within the errors hence the 
values quoted must be better considered as lower bounds. 



p 


6 


6 


6 


0.23142 


1.870(3) 


1.027(7) 


0.727(13) (?) 




(r = 3) 


(r = 3) 


(r = 3) 


0.2275 


2.593(4) 


1.429(8) 


1.016(13) (?) 




(r = 4) 


(r = 4) 


(r = 4) 



Table 6: Correlation lengths extracted with the variational approach in the 
3d 4> 4 model. 



P 


6 


6 


6 


0.405 


1.1130(9) 


0.609(2) 


0.464(5) 




(r = 2) 


(r = 2) 


(r = 2) 


0.385 


2.182(2) 


1.182(5) 


0.872(9) 




(r = 3) 


(r = 3) 


(r = 3) 


0.3798 


3.463(4) 


1.827(8) (?) 


1.246(19) (?) 




(r = 4) 


(r = 5) ' 


(r = 4) 



constructed the ratio of these masses with the lowest one. 

These adimensional ratios should approach a constant as (3 — > j3 c and be 
universal, namely they sould have the same value in the Ising and 4 models. 
This is clearly confirmed in Tab. 7 where we have reported these ratios for the 
two models together. A naive fit of these data (neglecting possible corrections 
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Table 7: Mass ratios for Ising and <p A models. 



Model 


P 


6 


6/6 


6/6 


4 


0.405 


1.1130(9) 


1.828(7) 


2.40(3) 


Ising 


0.23142 


1.870(3) 


1.821(15) 


2.57(5) (?) 


4 


0.385 


2.182(2) 


1.846(10) 


2.50(1) 


Ising 


0.2275 


2.593(4) 


1.815(13) 


2.55(4) (?) 


4 


0.3798 


3.463(4) 


1.895(10) (?) 


2.78(5) (?) 



to scaling and neglecting the values denoted with a question mark which must 
be considered only as upper bounds) gives the two results 

^ = 1.83(3) , ^ = 2.45(10) . (48) 

mi mi 

It is important at this point to discuss the relationship of these results with 
the perturbative expansion in the <p A theory discussed in Sec. 3.4. The first 
excitation that we find above the lowest one is below the pair production 
threshold: therefore it cannot be a perturbative effect related to cuts in the 
Fourier transform. Hence it is a new non-perturbative state which exists in 
the spectrum of the two models and cannot be seen within the framework of 
a perturbative analysis. 

On the contrary the second mass state could well be related to the cut. 
As a matter of fact, if we ignore the power tail in the functional expression of 
the cut and try to mimic it with a simple exponential we would exactly find 
a fictitious state with a mass of approximatively 2.4 times the lowest one M. 
This state could well be identified with the third mass that we measure in 
our analysis. 

7 Duality and the Glueball spectrum 

As already noticed in Sec. 2.1 the correlation length of the Ising model is re- 
lated by duality to the (inverse of the) mass of the + glueball of the 7L 2 gauge 
model. It is natural to conjecture that the new states that we have found 
in the spectrum of the Ising model are related by duality to the remaining 
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Fig. 5: Comparison between the spectrum obtained with the varia- 
tional method in the Ising model at (3 = 0.2275 and the values of the 
masses in the + family of the 3D TL<i gauge model at the dual coupling 
(3 = 0.74883. Diamonds, squares and triangles denote £i,£2 and £3 re- 
spectively, while pluses, crosses and stars denotes the + , +/ and + " 
glueballs. 

states of the + family. This conjecture is strongly supported by our data. 
In particular it is impressive to compare the values of £ e ff obtained in this 
paper as a function of the separation r with the corresponding observables 
for the glueball states + , +/ , + " obtained in Ref. |§. This comparison is 
shown in Fig. 5 for the value (3 = 0.2275. While for the first two values of 
r (which are dominated by the lattice artifacts) we find large discrepancies, 
for t > £ the two sets of data remarkably agree within the errors. This iden- 
tification is further supported by the values of the mass ratios extrapolated 
to the continuum limit, which again agree within the errors (see Tab. 8). 



25 



Table 8: Mass ratios for the Ising model and 7L 2 gauge model. The results 
for the gauge model are taken from Ref. ||/ 



Mass ratio 


Ising model 


7L 2 gauge model 


6/6 
6/6 


1.83(3) 
2.45(10) 


1.88(2) 
2.59(4) 



8 The spin-spin correlator revisited. 

Once the spectrum has been understood we can again address the behavior 
of the spin-spin correlator. In order to clarify our analysis we devote the 
following two sections 8.1 and 8.2 to a detailed discussion of G(r) in the 
particular case of the Ising model at /3 — 0.2275 where both data from 
Ref. |2|, H and from the present variational analysis exist. To allow the reader 
to reproduce our analysis we report in Tab. 9 the values of the connected 
spin-spin correlator which were obtained in Ref. (to which we refer for 
information on the parameters of the simulations and the algorithm used). 

8.1 Overlap amplitudes 

In Sec 4.2 we constructed an estimator R e ff for the universal ratio R. The 
slow convergence of R e ff to its asymptotic value, shown in Fig. 3, was ex- 
plained as originating from higher mass states in the spectrum and cut ef- 
fects. Now that the low-lying part of the spectrum has been determined by 
the variational method, the preasymptotic behavior of the correlation func- 
tion G(t) is under control. This information allows us to construct a new 
estimator of R with better convergence properties. Moreover, the presence 
of non-perturbative states in the spectrum makes it natural to define and 
study new universal quantities related to the overlap of the spin operator on 
these states, that are the exact analogs of R for the new states. 

As discussed in the previous section, our results for the spectrum show 
a non-perturbative state with a mass of about 1.8 times the fundamental 
one, and a third state that could be another non-perturbative state or a 
perturbative interaction effect (cut), or a superposition of the two. These 
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Table 9: The connected correlator G(t) at (3 = 0.2275. The data are taken 
from Kef. ||/. In the first column we report the distance r in units of the 
lattice spacing and in the second column the value of G{r) x 10 5 . 



T 


G(t) x 10 5 





152.30447(5) 


1 


98.57691(4) 


2 


65.20074(4) 


3 


43.57136(4) 


4 


29.29125(3) 


5 


19.76683(3) 


6 


13.37177(3) 


7 


9.06161(3) 


8 


6.14757(3) 


9 


4.17518(3) 


10 


2.84039(3) 


11 


1.93695(3) 


12 


1.32855(3) 



two possibilities correspond to two different ansatze for the functional form 
of G(t): The first one (that we shall call in the following the "three mass 
ansatz") is: 

G(t) = ciexp(-<r/fi) + c 2 exp(-r/£ 2 ) + c 3 exp(-r/£ 3 ) . (49) 

This choice does not make use of any perturbative information on the theory, 
and is based on the only assumption that higher states, which are certainly 
present, are beyond our resolution, i.e. their contribution is of the same 
order of magnitude as our statistical errors and hence cannot be taken into 
account. 

The second possibility (that we shall call in the following the "two mass 
plus cut ansatz" ) is based on the assumption that the third state we see in the 
spectrum is not a new state but the effect of the cut, that is a perturbative 
effect due to the self interaction of the fundamental state. 

G(r) = ci(exp(-T/6) + /cut(r/£i)) + ca exp(-r/6) • (50) 
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and f cut given by Ref. [§] 




oo 



e 



fcut[T) 



djJL 



(51) 



Here we are approximating the cut contribution with its one loop expres- 
sion. Again, we assume that higher states in the spectrum are beyond our 
resolution. 

We shall show below that the results obtained with these two choices 
essentially coincide. Thus, with the data at our disposal, we cannot select 
between the two scenarios, but at the same time we are sure that this sys- 
tematic uncertainty does not affect our results for £1, c\, £2 and Ci. 

With the constants q obtained in this way we may construct estimators 
for the universal quantities: we define 



with % — 1, 2, 3 or i = 1, 2 according to the ansatz. The functions q (t) are 
determined from the Monte Carlo value of G(r) for three or two nearby values 
of t and from the values of the masses as determined from the variational 
method. is our new, improved estimator of the universal ratio R, while 
R e 2 ff and R e / f are estimators for the new universal ratios associated with 
the new states. Let us discuss in detail the results obtained with the two 
ansatze. 

8.1.1 The three mass ansatz 

The results may be summarized as follows: 

• The asymptotic estimate of R is consistent with the one obtained in 
Sec. 4.2 with the unimproved estimator R e ff, but R\^ reaches its 
asymptotic value already for r ~ £ (see Fig. 6). 

• Despite the fact that the presence of the higher masses and the related 
new degrees of freedom slightly increases the systematic uncertainty in 
the function i? x (r) the error on the estimate of R is slightly smaller 
than that quoted in Tab. 2. This is due to the fact that we can extract 




(52) 
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our asymptotic estimate already at r ~ £, where the statistical un- 
certainties on G(t) are smaller. Since the major source of uncertainty 
comes from the error in the estimate of £1 and gives a contribution which 
increases linearly with r this improvement compensates the above un- 
certainty. Our final result is R = 2.055(15) (to be compared with the 
value R = 2.068(18) obtained with the unimproved estimator). 

• The function R\ (t) shows a stable behavior in the range £ < r < 2£ 
(see Fig. 7) and a reliable estimate of R2 can be extracted. The result 
is affected by errors (which are mainly due to the uncertainty in £1) 
larger than those which affect R. Our final result is R 2 = 0.45(8) . 

• The function c^{r) never stabilizes and it is impossible to extract a 
reliable asymptotic estimate for R 3 . 

8.1.2 The two mass plus cut ansatz. 

Again using as input parameters the values of £1 and £ 2 extracted from the 
variational analysis we may construct from pairs of nearby values of G(r) 
the functions R\ (r) and R% (r). We find the following results: 

• The behavior of (t) is essentially unaffected by the change of ansatz 
(see Fig. 6). Our best estimate for R x does not change. 

• The function R^\r) is more affected by the change, but in the region 
£ < r < 2£ where the data reach a stable plateau the new estimate 
agrees within the errors with those obtained with three masses (see 
Fig. 7). Our estimate in this case is i? 2 = 0.55(5) . 

Thus we may conclude that the interpretation of the third mass as a cut is 
compatible with the data and suggests a slightly higher value for R 2 . 

Trying to take into account the systematic error involved in the choice of 
the two ansatz we give as our final result: 

R = 2.055(15) #2 = 0.50(10) (53) 

which is a suitable combination of the two estimates. 
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2.35 - 
2.3 - 
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eff 




Fig. 6: The improved estimator R± (r) derived from the three mass 
ansatz (pluses) and the two mass plus cut ansatz (squares). Also shown 
for comparison is the unimproved estimator R e ff(r) defined in Sec. 4 
(diamonds). The data are at (3 = 0.2275 and the distance r is measured 
in units of £. 
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Fig. 7: Values of R2 obtained assuming the three mass ansatz (dia- 
monds) and the two mass plus cut ansatz (pluses) for (3 = 0.2275. The 
distance r is measured in units of £. 
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8.2 The behavior of £e//( r )- 

As a final test of our results, let us now turn our attention to the function 
Ceff( T ) an d test if, from the knowledge of the spectrum and the overlap 
constants, we are able to reproduce the observed behavior of £ e ff( T )- Since 
the two ansatze discussed above give essentially equivalent results we have 
chosen to do this test with the two mass plus cut ansatz only. We have 
constructed our best estimate of the correlator, inserting into the equation: 

G th (r) = c 1 (eM-r/Ci)+faat(r/^))+c 2 exp(-r/6) • (54) 



our best estimates for c±, c 2 , £1 and £ 2 reported in Eq. fl5"3|) and Tab. 5 respec- 
tively. Then we constructed: 

^ (r) = ]n(G th (r+l))-]n(G th (r)Y ^ 

The result is compared with the observed Ce//( T ) m Fig- 8- All the data for 
T > £ agree within the error bars. For comparison, we show also the purely 
perturbative prediction. 



8.3 Continuum limit. 

In Sec. 4.2 we have shown that the ratio R has a good scaling behavior 
and we extracted an estimate of its value in the continuum limit. It is very 
important to test if also the ratio R 2 has good scaling properties. To this end 
we have studied its value for two other values of (3: [3 = 0.2260, (3 = 0.2240, 
using again the data of Ref. [@. Let us briefly motivate this choice and the 
procedure we used: 

• We chose these two samples because for values of (3 higher than (3 = 
0.2275 the value of £2 becomes so small that it is impossible to observe 
its overlap with the techniques discussed above. 

• In the analysis we need the values of £1, £2 and £3. which we have not 
evaluated explicitly in this paper in the two cases (3 = 0.2260, (3 = 
0.2240. However, as shown in Sec. 7, duality allows us to identify £1, £2 
and £3 with the masses of the first three states in the + family of the 
glueball spectrum of the 7L 2 gauge model, which we evaluated, precisely 
for these values of /3, in Ref. ||. 
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0.5 1 1.5 2 2.5 3 3.5 

Fig. 8: Data for £ e //(r) at (3 = 0.2275 (taken from Ref The solid 
line corresponds to the function ^ t h dehned in Sec. 8.2 The dotted line 
corresponds to the function £ pert dehned in Sec. 3.4. The distance r is 
measured in units of £. 
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We performed the same analysis discussed in the previous section and found 
analogous results, that is: 

• The estimates of the ratio R obtained with the unimproved and im- 
proved estimators are consistent; the improved estimator allows one to 
reach the asymptotic value at much shorter distances. 

• For each (3, the values for R 2 obtained assuming the two extreme situ- 
ations: the three mass and the two mass plus cut ansatz agree within 
the errors in the range £ < r < 2£. 

• In this same range, the estimate of R 2 show a perfect scaling behavior as 
a function of (3. All the values obtained for the three different samples 
and the two possible ansatz agree within the errors and fall into the 
window R 2 = 0.5 ±0.1 (see Fig. 9). 

8.4 Universality. 

Another important issue is to see if we find in the 4 model the same values 
for R and R 2 that we obtained in the Ising case. To answer this question we 
analyzed the data for the two samples nearest the critical point, along the 
line discussed above. For r ~ £, choosing the three mass ansatz we obtained 
the results reported in Tab. 10. For the value of (3 nearest the critical point 
we find a perfect agreement with the Ising results. 

Table 10: R\ and R 2 in the 4 model, using to the three mass ansatz. 



p 


Ri 


R 2 


0.385 
0.3798 


1.981(30) 
2.047(20) 


0.75(40) 
0.59(20) 



9 Conclusions 

In this paper we have developed a new variational method to study the spec- 
trum of statistical models. We applied it to the three dimensional Ising and 
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Fig. 9: Values of R2 obtained assuming the three mass ansatz for 
(3 = 0.2275 (diamonds), (3 = 0.2260 (squares) and (3 = 0.2240 (triangles), 
compared with those obtained with the two mass plus cut ansatz (pluses 
for (3 = 0.2275, crosses for = 0.2260 and stars for (3 = 0.2240 . The 
distance t is measured in units of £. All the data are taken from Ref. ||/. 
The distance r is normalized, for each (3, to the asymptotic value To 
avoid confusion only the error bars for the (3 = 0.2275 data are reported, 
the other error bars are of the same size. The dashed line corresponds 
to the asymptotic value R2 = 0.5(1) quoted in Sec. 8.1. 
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4 models and succeeded in detecting two states beyond the lowest mass 
excitation. For these new states we could give a rather precise estimate of 
the mass and, for the lowest of them, also of the overlap constant. By com- 
paring our results with those published in Refs. 0, || we reach the following 
conclusions: 

• The state denoted in the paper as £ 2 , which is the first one above the 
lowest mass, has a truly non-perturbative nature. 

• The next state, denoted in the paper as £3, could again be a non- 
perturbative excitation, but the data are also compatible with a per- 
turbative origin: it could be the signature of the infinite set of states 
above the threshold for pair production of the lowest mass state. 

• The non-perturbative states offer a natural explanation for the dis- 
agreement between perturbative predictions and Montecarlo results for 
the universal quantities i/^nd and R. 

• We have been able to extract estimates for the adimensional ratios: 
£ 2 /£i, ^3/^1, R, R2 which show a good scaling behavior as functions of 

P- ' 

• Universality holds for the entire spectrum of the model. The states 
that we observe in the 3d Ising model appear with the same masses 
and overlap constants also in the </> 4 model. 

• Duality holds for the entire spectrum of the model. The states that we 
observe in the Ising spin model are related by duality to the glueballs 
of the + family of the 3D 7L 2 gauge model. 
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Appendix: Perturbative calculations at two 
loop order 



We consider the 3D Euclidean field theory defined by the action 



S= (fx 



td^t + ± (0 2 - v 2 ) 

which we want to treat perturbatively around the stable solution 

(f> — v . 

Therefore we define the fluctuation field (p: 

4> = v + ip 

in terms of which the action is 

S = So + Sj , 



2\ 2 



where 



So 
Si 
m 2 



1 
I 

gv 
3 



(fx 
d 3 x 

2 



-d^d^ip + —cp 
rn^ 3 + 9_ 4 



(56) 

(57) 
(58) 
(59) 

(60) 
(61) 
(62) 



From this expression we can read the Feynman rules in momentum space: 
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We will need the one and two point correlation functions of the tp field at 
two loop order; the diagrammatic expansions are: 
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and 
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All of the relevant integrals have been calculated in dimensional regulariza- 
tion by Rajantie in Ref. [T5|], where one can also find the explicit expression 
of the one particle irreducible part of the two point function. We obtain for 
d = 3 - 2e 
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(63) 



where /i is an arbitrary mass scale and u is the dimensionless bare coupling 
u = g/m; the two point connected function gives 
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with 



p +m H m" 

47T 



3m p 

1 arctan 

2p 2m 



38 



+ 



+ 



+ 



u 



-m 



1 1 fx 2 Ap 2 + 25m 2 m p 

, j lo^ ctrct cLii 

12e 6 9m 2 4(p 2 + 4m 2 ) p 3m 



16tt 2 

2p 4 + 35m V + 54m 4 
24p 2 (p 2 + 4m 2 ) 

,2 



log 1 + 



P 



9m 2 



+ 



p 2 + 10m 2 
4(p 2 + 4m 2 ) 

P 



3 m 2 p 3 , o m p 3 
-—r- arctan 1 — log 3 — arctan 1 A 

4 2m 8 p 2m 16 Vm 



9 -B (L 

2 Vm 



(65) 



where p = \fp 2 and 
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B(x) = 
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These correlation functions are used to fix three renormalization constants 
in order to compare the perturbative results with lattice data. We use the 
scheme introduced in Ref. |iT6fl , and define 
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G- l {p)=Ztf K + P 2 + 0(p 4 )) 

v R = Z R 1/2 (4>) , 



so that vr is to be identified with the magnetization and m R with l/^nd- 
Moreover it is useful to define 



ur = 3^- , 



corresponding to the scaling quantity: 

3X 



e 2nd M 2 ' 

where x is the susceptibility and M is the magnetization. 
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^From Eqs. (^,^) we find 
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and 



u R = u 



1 1 j yu 2 40957\ 
32 4tt ' 16vr 2 V 247 ~ 12 ° g 9^ + 55296 / 



31 u u A 
1 - — — + 



(73) 



The physical mass m v h is defined as the location of the pole of G(p) on 
the imaginary axis, that is by imposing 
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We find 
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so that, putting together Eqs. (0>[OIEi) we 
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that is Eq. fl22|). 

Finally, we are interested in the perturbative evaluation of the amplitude 
ratio R defined in Sec. 4, which is given by 
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The result is 
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that is Eq. 
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